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V- 

Abstract. We present a systematic investigation of rotation curves (RCs) of fully hydrodynamically simulated 
galaxies, including cooling, star formation with associated feedback and galactic winds. Applying two commonly 
used fitting formulae to characterize the RCs, we investigate systematic effects on the shape of RCs both by obser- 
vational constraints and internal properties of the galaxies. We mainly focus on effects that occur in measurements 
of intermediate and high redshift galaxies. We find that RC parameters are affected by the observational setup, 
like slit misalignment or the spatial resolution and also depend on the evolution of a galaxy. Therefore, a direct 
comparison of quantities derived from measured RCs with predictions of semi-analytic models is difficult. The 
virial velocity V c , which is usually calculated and used by semi-analytic models can differ significantly from fit 
parameters like V max or V op t inferred from RCs. We find that V c is usually lower than typical characteristic 
velocities derived from RCs. V max alone is in general not a robust estimator for the virial mass. 

Key words. Galaxies: kinematics and dynamics - Galaxies: spiral - Galaxies: structure 



1. Introduction 

Spatially resolved rotation curves (RCs) are a fundamen- 
tal tool to study the internal kinematics and the distri- 
bution of mass in spiral galaxies. The first discrepancy 
between the differential kepler-type of rotation and the ro- 
tation of galaxies was detected by Babcock (1939). These 
first indirect measurements of non-luminous matter intro- 
duced the concept of dark matter (DM) into astrophysics. 
An overview is given by Sofue and Rubin (2001) and ref- 
erences therein. 

An important application of RCs lies within a correlation 
of the luminosity and the maximum rotational velocity of 
spirals found by Tully & Fisher (1977). The physical origin 
of the slope and the scatter of the TFR is still subject to 
debate. Different theoretical approaches exist, which dif- 
fer mainly in the predictions of the redshift evolution of 
the TFR. Therefore, Ziegler et al. (2002) and B6hm et al. 
(2004) used a sample of field galaxies in the FORS Deep 
Field to study the TFR at intermediate redshift. They find 
a significant change of slope in comparison to local sam- 
ples, mainly caused by small, star forming distant galaxies. 
However, the measurement of rotational velocities is more 
complicated in the case of distant, apparently small spi- 
rals. 

In this work we investigate how parameters from models, 
describing the shape of RCs, are influenced by observa- 
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tional constraints and by internal properties of galaxies. 
Only techniques that take these systematics into account, 
as e.g. the method presented in Bohm et al. (2004), can 
get robust results for V max . 

In recent years, fully N-body/hydrodynamic simulations 
of spiral galaxies became an important tool to understand 
the formation and evolution of spiral galaxies (e.g. Mihos 
& Hernquist 1994; Springel & Hernquist 2002). In these 
simulations cooling, stellar feedback and galactic winds 
are taken into account to model galaxies in a physically 
motivated way. Here, we extract RCs from model galaxies 
simulated with GADGET2 (Springel 2005). We investi- 
gate for intermediate and high redshift galaxies the influ- 
ence of large relative slit widths, inclinations and slit mis- 
alignments on the determination of fitting parameters like 
V op t or V max - These parameters are commonly used (e.g. 
Courteau 1997, Yegerova & Salucci 2004) as a measure for 
the 'peak' circular velocity, e.g. to determine Tully-Fisher 
relations. 



2. Simulations 

Recently Kapferer et al. (2005) studied the influence of 
galaxy-interactions on the strength and evolution of the 
star formation rate of the interacting system. In this work 
we investigate the isolated model galaxies, presented in 
Kapferer et al. (2005). The initial conditions (hereafter 
ICs) of the model galaxies were built according to Springel 
et al. (2004), which is based on the work of Mo et al. 
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Table 1. Properties of the initial conditions of the model 
galaxies 



Properties 


Galaxy A 


Galaxy B 


circular velocity V C L 


160 


80 


disk mass fraction 2 


0.05 


0.05 


gas content in the disk 3 


0.25 


0.25 


disk thickness 4 


0.02 


0.02 


total mass [Mq] 


1.33xl0 12 hr 1 


1.67xlO n h' 1 


disk scale length [kpc] 


4.51k- 1 


2.25k- 1 



circular velocity at r2oo in km/s 
... fraction of disk particles (stars/gas) in units of halo mass 
3 ... relative content of gas in the disk 
4 ... thickness of the disk in units of radial scale length 



(1998). The two model galaxies in this work are chosen 
such that they represent a Milky Way type and a small 
spiral galaxy without any bulge component. In Table^the 
properties of the ICs of the model galaxies are listed. The 
combined N-body/SPH simulation calculates then 5 Gyr 
of isolated evolution. For every time step we know the 
velocity of each particle and can hence extract realistic 
rotation curves. 

2.1. Rotation Curve Extraction 

In order to extract the rotation curves (RCs) of our sim- 
ulated model galaxies, we define a slit with a width d, 
see Fig. ^ In addition we allow for a misalignment an- 
gle S to simulate rotations of the slit with respect to the 
major axis of the system. Such rotations sometimes occur 
in observations using multi-object spectroscopy. In Fig. ^ 
the different parameters for the slit are shown. The slit 
width d and the misalignment 5 of the slit with respect to 
the major axis define the virtual slit. We extract the RCs 
from the velocity field of the gas in the following way. The 
velocity as a function of radius is determined by averag- 
ing over all line-of-sight velocities of gas particles in thin 
bins, with a bin size Ar along the slit and a side length 
d perpendicular to the spatial axis. In Fig. [21 the veloc- 
ity field taken for the RC extraction is shown (a). In Fig. 
(b) the spatial sampling along the slit is sketched. An RC 
with a resolution of 0.1 kpc was extracted and used as a 
reference. 

In order to determine systematic effects of large slit widths 
d (relative to the galaxy size), which occur in the case of 
observations of distant (z«0.5) galaxies (e.g. Bohm et al. 
2004) we vary d in a range of several kpc. Not only the 
relative slit width varies in observations of distant disk 
galaxies but also the sampling of the velocity field along 
the spatial axis (i.e. the spatial resolution). In order to 
simulate this finite spatial resolution we bin the reference 
RC with different bin sizes. In Fig. [3] we show the extrac- 
tion of different RCs, corresponding to different spatial 
resolutions along the slit. 




Fig. 1. Image of model galaxy A and a virtual slit for 
extracting a rotation curve. The slit width d and the slit 
misalignment angle 5 are indicated. 





ISHH||HH 







Fig. 2. Sketch of our procedure to extract rotation curves 
of the model galaxies, (a) the line-of-sight velocity field 
of a model galaxy is indicated as it would be observed 
by a virtual observer, (b) the sampling along the slit is 
highlighted (Ar). 



3. Results 

3.1. Rotation curves as a function of evolution 

As a first step we investigate the RCs of the ICs and 
compare them to the RCs of the fully hydrodynamically 
treated galaxies after 5 Gyr of evolution. Note that the 
ICs are based on an analytic model introduced by Mo et 
al. (1998). The evolution starting from these ICs is de- 
termined by the influence of the dynamics, the star for- 
mation with feedback and stellar winds of the system. 
Therefore the RCs and other internal properties of the 
evolved galaxies do differ from those of the ICs. The gen- 
eral evolutionary trend of the RCs is presented in Fig. 
The rotation curves were obtained by setting a slit width 
of 4 kpc (galaxy A) and 1 kpc (galaxy B) without any slit 
misalignment <5. The galaxies were inclined with an incli- 
nation angle of i=80°. A spatial resolution of 0.1 kpc was 
adopted to extract ideal RCs. It is clearly visible, that the 
rotational velocities get lower for the evolved galaxies. The 
decrease of the overall angular momentum of gas particles 
in the disk can be explained by mass ejection due to galac- 
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Fig. 3. An example for a 5 kpc binning of the ideal RC. 



tic winds and the rearrangement of the gas in the disk due 
to the fully hydrodynamic treatment as the disk rotates. 
As we want to investigate the dependencies on different 
observational constraints, e.g. galaxy alignment with re- 
spect to the spectroscopic slit or slit misalignment, we use 
hereafter models to describe the shape of the RC. Used in 
observational work are physically motivated fitting func- 
tions, like the universal rotation curve (URC) eq. 2 (Persic 
et al. 1996), or purely phenomcnological fitting formulae 
like eq. 1 (Courteau 1997). Although we are aware that 
eq. 1 cannot reproduce the many observed declining RCs, 
it is suitable for our model galaxies, which do not decline 
in the outer parts. The RCs shown in Fig. 0] are best fits 
to the measured rotational velocities of the model galaxies 
using the fitting formula of Courteau (1997). This function 
is defined as 



V rot (r) 



(r a + rg)a 



[km/s] 



(1) 



where r is the galactocentric distance and ro and a are 
free fitting parameters. Clearly, with the general decrease 
of the rotational velocities also V max gets lower. The phys- 
ically motivated URC (Persic et al. 1996) is a superposi- 
tion of the velocity field of the DM halo and the disk. The 
URC can be expressed as follows 



Uurc I — ) = V(r opt )\ f 0.72 + 0.44 log(^) 



opt 



1 97x 122 
(^ + 0.78^)^3 +L6exp [-°- 4 



L 



x 2 + 1.5 : 



[km/s], (2) 



where x — r/r op t, V op t is the rotational velocity at r op t 
(the reference scale, which for an exponential disc is 
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Fig. 4. Best fits to the measured RCs of model galaxies A 
and B. The fits are done for the ICs and for the evolved 
(5 Gyr) systems. 
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URC (Persic, Salucci & Stel 19 
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Fig. 5. Comparison of the two profiles (Courteau, URC) 
for model galaxy A. Both profiles represent the velocity 
field very well. 



3.2i?d) and L the absolute blue luminosity of the galaxy. 
In Fig. [5] the URC fit to our data is shown together with 
the Courteau (1997) fit. Both functions represent the RCs 
of the model galaxies very well. Although we are aware 
that the Courteau fitting function has no physical justifi- 
cation, we use both functions to fit our data, in order to 
investigate the shape of the RC as a function of the obser- 
vational setup, i.e. galaxy alignment and slit properties. 
The fit parameters for eq. 1 and 2 with 95% confidence 
level errors are listed in Table Eland 3, respectively. Note 
that no restrictions were set on the fit parameters (V max , 
ro and a for the Courteau function and V op t for the URC). 
The blue band luminosity L b was estimated from the stel- 
lar mass assuming a stellar mass to light ratio of 1.2 (mean 
from different star formation histories in the redshift range 
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Table 2. Fitting parameters for eq. 1 (Courteau 1997) 
for RCs of the ICs and the fully hydrodynamically treated 
galaxies (evolution 5 Gyr). 



Fitting 


Parameter 


Galaxy A 


Galaxy B 


I 
I 
I 


Vmax [klTl/s] 

ro [kpc] 
a 


205+a 

Qi+0.14 
°- a± -0.14 

1 62+ 011 


1 OQ+0.07 
i.oo_ og 

1 89+ ' 51 


II 
II 
II 


V max [km/s] 

ro [kpc] 
a 


169.67+457 
4 6+ - 4 

^■"-0.47 
o oq+1.24 
o.o»_ 1 23 


31, °-2.36 

2 05+° 17 

Z.UJ_ 25 

2 14+ - 41 


I ... initial conditions 



II ... evolved galaxies (t=5 Gyr) 



Table 3. Parameter V op t for eq. 2 (Persic et al. 1996) for 
RCs of the ICs and the fully hydrodynamically treated 
galaxies (evolution 5 Gyr). 



Fitting 


Parameter 


Galaxy A 


Galaxy B 


I 
I 


Vopt [km/s] 
L/L, 


186.5+^ 
0.61 


80.3+^ 
0.075 


II 
II 


Vopt [km/s] 
L/L* 


183.2+g 1 
0.65 


91-1-2:2 

0.08 


I ... initial conditions 



II ... evolved galaxies (t=5 Gyr) 



z = 0.5 — 1.4, Dickinson et al. 2003). The L* luminosity 
at z = 0.5 was adopted as M% = —21.3 (Gabasch et al. 
2004) this corresponds to log(L* B ) = 10.72. We find val- 
ues for L/L* on the range of 0.61-0.65 for model galaxy A 
and 0.075-0.08 for model galaxy B. An evolution of L/L*, 
due to new forming stars in the galaxies, can be seen. For 
the investigations of the RC as a function of observational 
bias, we choose L/L* as an additional free parameter, to 
get better representations of the RCs. 
Most present models of galaxy formation and evolution 
rely on the work of Mo et al. (1998). Indeed, this model can 
reproduce correctly the general shape of rotation curves, 
but galaxy evolution can alter the RC. Note especially 
the differences between V c and V max / opt from tables |21 
and 3, respectively. A comparison of semi-analytic models 
and observations is generally complicated by the fact that 
V max /opt is determined differently from V c . In the case 
of mass reconstruction via RCs the superposition of the 
velocity field of the halo and the disk in the URC ansatz 
would be the adequate approach. 

3.2. Rotation curves as a function of the slit width for 
model galaxy A 

In order to study systematic effects of large relative slit 
widths, as they appear in observations of galaxies at in- 
termediate and high redshift, we extract the RC of model 
galaxy A for several slit widths. Of course slit widths for 
local spiral systems are orders of magnitudes smaller. But 
for observations of galaxies in the redshift range 0.5 to 
1, as carried out in a project related to the present work 



250 r 




-250 1 1 1 1 1 1 1 ' 

-20 -15 -10 -5 5 10 15 20 

radius [kpc] 

Fig. 6. RCs of model galaxy A ICs, for two different slit 
widths (1 kpc and 10 kpc). The error bars indicated in the 
figure are the standard deviations of the mean velocity in 
each 0.1 kpc bin, see Fig. |3] (image b). 



(Ziegler et al., 2003), typical slit widths become compara- 
ble to the disk scale length Rd. It is important to note, that 
such large slit widths result in an integration of the veloc- 
ity field perpendicular to the spatial axis (slit direction). 
This effect is the optical equivalent to 'beam smearing' 
in radio observations. To investigate this effect we mea- 
sure RCs for slit widths ranging from 1 kpc to 10 kpc, 
which corresponds for our model galaxy A to 0.2ft, -1 Rd 
to 2.2fr~ 1 Rd. In Fig. El and [7| we show the RCs for model 
galaxy A (IC and evolved, respectively) for two very dif- 
ferent slit widths (1 kpc and 10 kpc). 
If the slit width is 1 kpc the scatter around the mean ve- 
locity in each bin is very small in comparison to the 10 kpc 
slit. This can be explained in terms of velocity distribu- 
tions in a bin. As the slit width is increasing more gas par- 
ticles can contribute to the measured mean velocity in a 
bin. In other words, the mean velocity is a superposition of 
velocity components from different regions of the galaxy, 
mainly due to the line of sight velocity distribution. In Fig. 
0the same quantity is shown as in Fig. El but after 5 Gyr of 
evolution of model galaxy A. Again the same behaviour in 
the scatter and mean velocity is present, but the overall ve- 
locity field shows more structure. This is a consequence of 
the fully hydrodynamic treatment of the galaxy. Note that 
ICs (as analytic models) do not include prescriptions for 
spiral arms, which are present in observed galaxies. Only 
the N-body/SPH simulations can reproduce this feature. 
Therefore the measured ideal RC for the evolved galax- 
ies shows local fluctuations, connected to e.g. spiral arms. 
This fact is well known from observations, where fluctua- 
tions of a few tens of km/s are superposed on the smooth 
rotation curve of the galaxy due to spiral arms (see e.g. 
Sofue and Rubin, 2001). This is in good agreement with 
our model RCs (cf. Fig. El and Fig 01. 
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Fig. 7. RCs of model galaxy A after 5 Gyr of evolution, 
for two different slit widths (1 kpc and 10 kpc). The error 
bars indicated in the figure are the standard deviations of 
the mean velocity in each 0.1 kpc bin, see Fig. |3] (image 
b). 



As a next step we fit eq. 1 to the RCs, extracted from 
different slit widths, shown in Fig. |SJ We used galaxy 
model A after 5 Gyr of evolution. The galaxy was always 
'observed' with an inclination i = 80°, with an ideal spa- 
tial resolution of 0.1 kpc. As the slit width increases the 
fitted curves decrease. Again this can be explained by the 
averaging process. It is obvious that too wide slits (d>Rd) 
result in non-flat RCs, which should not be fitted by eq. 
1. Instead an observer would use here the URC. If we 
adopt our fitting procedures and extract V max and V op t 
for different slit widths, we obtain a dependence of V ma x 
and V op t on d, as shown in Fig. |5J A nearly linear de- 
crease of V max and V op t from d=l kpc to d=4.5 kpc can 
be seen. From this investigation we would recommend to 
apply only a maximum slit width in the order of Rd . 



3.3. Rotation curves as a function of inclination 

Galaxies are very rarely observable edge-on and therefore 
a correction of inclination effects on the RC is important. 
The intrinsic RC V mt (r) of a galaxy is most often cor- 
rected by the sine of the inclination angle, i.e. by the sim- 
ple geometric correction V lnt (r)=V obs (r)/sin(i) (edge-on 
galaxies are defined to have i=90°). We investigate the 
RC of model galaxy A by rotating the galaxy from i=20° 
to i=80° , which is a typical range of inclination angles ac- 
cessible in observations. The slit was centred on the gas 
disk with a slit width fixed to d=4 kpc. After each rotation 
we extract the RC, corrected with the simple expression 
V int (r)=V obs (r)/sin(i) before, and fit the data with eq. 1 
and 2. The dependence of V max and V op t on the inclina- 
tion angle i is shown in Fig. 1101 In the inclination range 
80° > i > 70° V max and V opt show a steep increase while 
for i < 65° V max and V op t stay roughly constant. The ex- 




- slit width 1 kpc 
slit width 3 kpc 
slit width 5 kpc 
slit width 7 kpc 



8 10 12 14 16 
radius [kpc] 

Fig. 8. Best fits to RCs for different slit widths. The un- 
derlying galaxy model is A after 5 Gyr of evolution. 



200 




3 4 

slit width [kpc] 

Fig. 9. Fitting parameters V max (blue) and V op t (red) as a 
function of slit width for model galaxy A with ideal spatial 
resolution. The dashed lines correspond to the bounds of 
the 95% confidence level. After an almost linear decrease 
the scatter increases for V max . 



planation for this behavior is shown in Fig.^] If a galaxy 
is observed nearly edge-on, the rotational velocity is an 
average of velocities from all radial distances along the 
line-of-sight. With decreasing inclination, more and more 
gas particles with lower line-of-sight velocity components 
move out of the slit, and are therefore not taken into ac- 
count in the averaging process. Thus, the mean velocity in 
each bin increases, which again leads to a larger V max and 
V op t- Below a certain inclination angle, depending on the 
slit width, most of the volume of the disk is not covered 
by the slit. In this volume most gas particles with a low 
line-of-sight velocity component are located. In our case 
70% of the volume is not covered by the slit for i < 65°. 
Note that this behavior is most significant for large rel- 
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30 40 50 60 70 80 

inclination angle (i) 

Fig. 10. Fitting parameter V max (blue) and V op t (red) 
as a function of the inclination angle. The dashed lines 
correspond to the bounds of the 95% confidence level of 
Knax- The underlying model is galaxy A, after 5 Gyr of 
evolution. A constant slit width of 4 kpc was used. 




Fig. 11. Sketch of the influence of the inclination on a 
spiral galaxy for a fixed slit width for model galaxy A. 



ative slit widths. In observations one tries to overcome 
the problem occurring at high inclination angles by spe- 
cial techniques, as e.g. the 'envelope-tracing' method (e.g. 
Sofuc and Rubin, 2001). 

As the inclination angle i decreases the errors for the fit 
become larger. Note that for i < 35° the errors are in the 
order of 10%. In Fig. El we show the fitting parameter 
V max as a function of the inclination angle i for different 
slit widths. If the inclination is below 70° and above 30° 
the slit width does not affect V max . The same behaviour 
was found for V op t. Only in the cases near edge on and 
face on, the slit width plays an important role in covering 
gas particles. The overall trend is the same as shown in 




175 - 



— slit width 2 kpc 
slit width 3 kpc 

170- slit width 4 kpc 

slit width 5 kpc 



165 L 



40 50 60 

inclination angle i 



Fig. 12. Fitting parameter V max as a function of the in- 
clination angle, for slit widths in the range of 2<d<5 kpc 
for model galaxy A. 



Fig. El where larger slit widths lead to lower V max and 
larger errors. 

Therefore, for large inclinations (i > 70° the slit width 
varies V max strongly (see Fig.lllland ll2l) . while for smaller 
inclinations the determination of V max does not strongly 
depend on the slit width. 

3.4. Rotation curves as a function of misalignment 

As shown in Fig. ^ the slit for measuring the RCs can 
have a misalignment S, with respect to the major axis of 
the projection of the galaxy on the sky. Especially multi- 
object spectroscopy has to deal with misaligned slits, 
therefore we investigate the effect of i5 on the fitting pa- 
rameters V max and Vopt- In Fig. ^1 we plot the result 
for different S for a fixed slit width of d=4 kpc for model 
galaxy A. Note that we applied the standard cos _1 (<5) cor- 
rection. Nevertheless the fitting parameters V max and V op t 
are not independent of 5. The corresponding error (scat- 
ter around the mean velocity in a bin), does not show 
any dependence on the misalignment angle i5, if it covers 
particles all along the slit. As the cos (5) correction, that 
we have applied, is only fully valid for two dimensional 
disks, without any thickness, V max and V op t show a de- 
pendence on <5, the slit misalignment. In the case of multi- 
object spectroscopy, where the misalignment can be much 
higher, more advanced corrections have to be applied. We 
have introduced one method in Bohm et al. (2004). 

3.5. Rotation curves as a function of binning 

To simulate different spatial resolution we bin our ideal 
RC, see Fig. |3 In Fig. ^] we show the dependency of the 
fitting parameters Kn ax /opt on different spatial resolutions, 
together with standard deviations. Obviously a poor spa- 
tial resolution leads to larger errors, but the value V max 
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Fig. 13. Fitting parameters V max (blue) and V op t (red) 
as a function of the misalignment angle, for a fixed slit 
width of d=4 kpc for model galaxy A. The dashed lines 
correspond to the bounds of the 95% confidence level for 
the fitting parameter V^ lax . 




2 2.5 3 
binning [kpc] 



Fig. 14. Fitting parameters V max and V op t as a function 
of the spatial resolution. Obviously the quality of an RC 
decreases for poorer spatial resolutions. 



shows no dependence. As V max represents the flat part of 
the RC (the asymptotic velocity), the binning does not 
vary V max , but of course a and rn. 

3.6. Estimators for the virial mass 

One major goal of investigating RCs is the possibility to 
determine the mass of the system including baryonic and 
non-baryonic components. The virial mass of a galaxy is 
cx to Vj? (e.g. Mo et al. 1998), where V c is the rotational 
velocity of the galaxy at the virial radius. Unfortunately, it 
is not possible to measure this velocity directly, therefore 
the virial mass of a galaxy is estimated by some suitable 



measure of the maximum circular velocity. If V max con- 
verges to V c at the virial radius, the RC fitting procedure 
with eq. 1 is a good estimator for the virial mass. However, 
as V max depends on observational constraints and of in- 
ternal properties of the galaxy, it seems to be no robust 
estimator for the total virial mass. Van den Bosch (2002) 
investigated the impact of cooling and feedback on disk 
galaxies using an analytical model. The author comes to 
similar conclusions, however, as he does not treat the gas 
self consistently (i.e. fully hydrodynamically) , he cannot 
derive a rotation curve and V max as we do. Instead in van 
den Bosch's work V max is defined as the maximum rota- 
tion velocity inside the radial extend probed by the cold 
gas (van den Bosch, 2002). The author finally concludes 
that eq. 3 seems to be a good fitting function to obtain 
the virial mass of a system. 



M vir = 2.54 x 10 10 Mfi 







I'd 
kpc 



(V max 
lOOkms- 



(3) 



Van den Bosch (2002) did not intend to give an exact es- 
timator of the virial mass by this fitting function. Instead, 
it seems that he wanted to give a rough estimate for the 
influence of cooling and feedback on the determination of 
the virial mass using V max . Therefore, it is clear that eq. 
3 is not exact. Nevertheless, as some observers use this 
function we also apply it to our simulated galaxies. 
As we treat the galaxies fully hydrodynamically, we are 
able to show the influence of evolution on the RCs. In 
Fig. \M the general decrease of V max is given. The main 
part, which decreases the angular momentum is the pres- 
ence of a galactic wind, which ejects matter from the disk 
into the surrounding halo. In most cases we find that V max 
derived from the fitting procedures does not represent V c 
and is also different from the maximum rotation veloc- 
ity present in the full velocity field. Most determinations 
of V max (different inclination, slit misalignments and slit 
widths) result in an overestimation of V c , if the standard 
corrections sin(i) and cos(<5) are applied. The virial mass of 
our model galaxies, obtained by extracting r2oo as limiting 
radius, yields to M2oo=l-13x 10 12 Mq for model galaxy A 
and M 2 oo = l-44x 10 11 Mq for model galaxy B. We used eq. 
2 to determine the virial mass and find a systematic under- 
estimation of w 50%. As van den Bosch (2002) states cor- 
rectly, matter can be ejected by galactic winds and there- 
fore reduce V max by increasing at the same time R^. The 
same behavior is present in our N-body/SPH simulations. 
However, we find that eq. 3 underestimates the virial mass 
of our model galaxies. The same underestimation of van 
den Bosch's (2002) estimation was stated by Conselice et 
al. (2005). However, as we have mentioned earlier, eq. 3 
is not thought to be exact. In fact, van den Bosch (2002) 
mentions that the error for an individual galaxy can still 
exceed a factor of 2. Thus it is no surprise, that the result 
is not correct for our model galaxy. 

A more detailed mass decomposition by applying the URC 
fitting would allow a deeper insight into the mass distri- 
bution of a galaxy. Nevertheless, the same problems from 
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the observational point of view would be inherent, like slit 
width, slit misalignment or galaxy orientation. A detailed 
investigation of the influence of the previous mentioned ef- 
fects on the mass decomposition by the RC will be investi- 
gated. Here we conclude that it is important for observers 
to investigate the environment of the measured galaxy. 
As galactic winds can strongly influence the internal kine- 
matics of the gas in the disk, the knowledge of the star 
formation rate gives constraints on the robustness of the 
determination of the virial mass from V max . An important 
issue in this context is the membership of the galaxy to a 
group, a galaxy cluster or the field. It is important to note 
that especially spiral galaxies in galaxy clusters can often 
interact with each other. As the interaction (Kapferer et 
al. 2005) increases the star formation rate significantly, 
merger-driven starbursts occur for a short (up to several 
100 million years) time and expel huge amounts of inter 
stellar matter into the surrounding halo. Therefore, these 
systems might have lower rotational velocities than iso- 
lated, low star forming galaxies while the general shape of 
the RC remains similar. 

3.7. The Tully-Fisher relation of spiral galaxies 

Semi-analytical galaxy formation models usually fix their 
free parameters such that the models match observed 
present-day luminosity functions (LF) or TFRs. Early 
models had problems predicting at the same time the cor- 
rect LF and the correct zero point of the TFR. Therefore, 
additional physical processes were introduced to fit both 
important statistical properties of a galaxy population at 
the same time. However, most semi-analytical models only 
have rough models to approximate the circular velocity of 
the halo V c , to derive a TFR for the underlying galaxy 
population. This in turn, can differ significantly from the 
maximum rotational velocity V max , derived from observed 
RCs, as discussed in the previous sections. Thus, differ- 
ences of the zero point of observed and modelled TFRs 
could, at least to some extend, origin from these discrep- 
ancies of V max and V c . We do not claim here, that the 
differences between V max and V c solve the problem of the 
zero point discrepancies of observed and simulated TFRs. 
We just point out that these two quantities are in gen- 
eral not equal and therefore it is problematic to use them 
equivalently. A major contribution to the difference be- 
tween V max and V c are galactic winds, as they are able to 
decrease angular momentum of a disk, by expelling a sig- 
nificant amount of matter into the surrounding halo. As 
we are now able to show, that for given modelled galaxies, 
V max as derived in observations, indeed mostly overesti- 
mates V c , we emphasise the importance of cooling and 
feedback processes. We plan to investigate the TFR of 
a semi-analytical model, taking this discrepancy into ac- 
count, in a forthcoming work. The effect of the overesti- 
mation of V max as V c is sketched in Fig. ^j] If V max was 
equal to V c we would expect that observed and simulated 
TFRs coincide in the right line. If V max is higher than 




1 toflO^J [knVs] 

Fig. 15. Sketch of a shift in the Tully-Fisher relation 
due to discrepancies between V c and V max . Only if the 
assumption V max =V c holds, Tully-Fisher relations from 
semi-analytical models and from observations are compa- 
rable. If V max >V c (result of our investigation) the line of 
the semi-analytical model would be shifted to the left with 
respect to the observed line. 



V c the TFR of the simulations would be located at the 
position of the left indicated line. 

It is important to stress, that the direct comparison 
of observation (V max ) and simulation (V c ) does invoke 
uncertainties originating from internal kinematics of the 
disk, like galactic winds. Similar concerns were empha- 
sized by van den Bosch (2000). Tully Fisher relations can 
also be constructed by using V op t ■ Again the same observa- 
tional constraints would introduce uncertainties, but the 
increase of the disk scale length due to galactic winds is 
compensated by measuring V op t at a given optical radius. 
Different disk scale lengths result in different optical radii 
and therefore the evolution of the gaseous disk is taken 
into account. 

4. Summary and conclusion 

In this work we use fully hydrodynamically modelled 
galaxies, including star formation, stellar feedback and 
galactic winds to study the internal kinematics of the gas 
in a spiral galaxy. We extract an RC from the line-of- 
sight velocities of the gas particles. Then we use a three- 
parameter fitting formula to describe the rotation curve 
(Courteau, 1997) and the universal rotation curve (URC, 
Persic et al. 1996). We find 

— that for evolved model galaxies the extracted and fit- 
ted RCs show a tendency to lower rotational veloci- 
ties, compared to the initial conditions. This can be 
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explained by galactic winds, which expel a certain 
amount of matter into the surrounding halo and there- 
fore decrease the total angular momentum. 

— We show that the variation of the slit width does influ- 
ence the quality of the RC and the values of V max as 
well as V pt- If the slit width is small, V max and V^, p t 
are systematically higher and show less scatter. If the 
slit width is large, the mean velocity is a superposition 
of particles along the line of sight, which introduces 
low velocity components and therefore larger scatter. 

— The dependence of V max and V op t on the inclination 
angle is nearly constant over the range of 20° < i < 
70°. In the range above i = 70° both show a strong 
dependence on i, which can be explained by rotating 
velocity components from the foreground and back- 
ground of the disk into the slit. It turned out that 
the sine correction leads to similar results in the range 
of 20° < i < 70°. In this range we do not encounter 
strong dependencies on the slit width. 

— The cos (8) correction for the slit misalignment is only 
fully valid for two dimensional disks, without any 
thickness. Thus, V max and V op t show an almost lin- 
ear dependence on S. In the case of multi-object spec- 
troscopy, where misalignments are inherent, more ad- 
vanced corrections have to be applied, as for example 
introduced by Bohm et al. (2004). 

— The spatial resolution does not influence V max and 
V pt strongly, but influences the quality of the RC. 

— We test the capability of V max as an estimator for the 
virial mass of the system and found a strong overesti- 
mation of the virial mass, by applying the virial theo- 
rem. By testing a more sophisticated relation, includ- 
ing results of semi-analytic models, introduced by van 
den Bosch (2002), we find an underestimation in the 
order of 50% of the virial mass. The explanation for 
the disagreement with van den Bosch lies in the fully 
N-body/hydrodynamic treatment in our simulations. 
Another point is the discrepancy of deriving V max , in 
our case from RCs from our model galaxies and in his 
case of the semi-analytical approach. 

— As Tully-Fisher relations are a common tool for testing 
models of galaxy evolution, any systematic differences 
between observations and theory play an important 
role. We show that V max and V opt usually differ from 
V c , which introduces a shift into the Tully-Fisher re- 
lation. 

The investigation of the influence of minor/major 
mergers and galaxy flybys is an important issue on this 
topic, which will be discussed in a forthcoming paper. 
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